Parametric coupling for superconducting qubits 
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We propose a scheme to couple two superconducting charge or flux qubits biased at their symmetry 
points with unequal energy splittings. Modulating the coupling constant between two qubits at the 
sum or difference of their two frequencies allows to bring them into resonance in the rotating frame. 
Switching on and off the modulation amounts to switching on and off the coupling which can be 
realized at nanosecond speed. We discuss various physical implementations of this idea, and find 
that our scheme can lead to rapid operation of a two-qubit gate. 
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The high degree of control which has been achieved 
on microfabricated two-level systems based on Joseph- 
son tunnel junctions [H [3- C3l has raised hope that they 
can form the basis for a quantum computer. Two ex- 
periments, representing the most advanced quantum op- 
erations performed in a solid-state environment up to 
now, have already demonstrated that superconducting 
qubits can be entangled [J, |5| • Both experiments imple- 
mented a fixed coupling between two qubits, mediated 
by a capacitor. The fixed-coupling strategy would be 
difficult to scale to a large number of qubits, and it is 
desirable to investigate more sophisticated schemes. Ide- 
ally, a good coupling scheme should allow fast 2-qubit 
operations, with constants of order 100MHz. It should 
be possible to switch it ON and OFF rapidly with a high 
ON/OFF ratio. It should also not introduce additional 
decoherence compared to single qubit operation. Charge 
and flux qubits can be biased at a symmetry point 0, l| 
where their coherence times are the longest because they 
are insensitive to first order to the main noise source, 
charge and flux-noise respectively. It is therefore advan- 
tageous to try to keep all such quantum bits biased at 
this symmetry point during experiments where two or 
more are coupled. In that case, the resonance frequency 
of each qubit is set at a fixed value determined by the 
specific values of its parameters and can not be tuned 
easily. The critical currents of Joscphson junctions are 
controlled with a typical precision of only 5%. The charge 
qubit energy splitting at the symmetry point depends lin- 
early on the junction parameters so that it can be pre- 
dicted with a similar precision. The flux-qubit energy 
splitting (called the gap and noted A) on the other hand 
depends exponentially on the junctions critical current 
9] and it is to be expected that two flux-qubits with 
nominally identical parameters have significantly differ- 
ent gaps 0- Therefore the problem we would like to 
address is the following : how can we operate a quan- 
tum gate between qubits biased at the optimal point and 
having unequal resonance frequencies ? 

We first discuss why the simplest fixed linear cou- 
pling scheme as was implemented in the two-qubit ex- 
periments 0, Q fails in that respect. Consider two flux 
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FIG. 1: (Color online) a) Two flux-qubits (shown coupled 
to their read-out SQUIDs and to their flux control line d 
(i = 1, 2)) coupled by a fixed mutual inductance M. b) Para- 
metric coupling scheme : the two flux-qubits are now coupled 
through a circuit that allows to modulate the coupling con- 
stant g through the control parameter A. 



qubits biased at their flux-noise insensitive point jq = tt 
(jq being the total phase drop across the three junc- 
tions), and inductively coupled as shown in figure 7]. 
The uncoupled energy states of each qubit are denoted 
|0j),|lj) (i = 1,2) and their minimum energy separation 
hAi = hiuJi. Throughout this article, we will suppose 
that Ai > A 2 . As shown before 0,0 , the system hamil- 
tonian can be written as H = H q \ + H q2 + Hj, with 
H qi = —(h/2)Ai<7 z i (i = 1,2) and Hi = hg o x io- x2 = 

l i go(< J t (7 2 + a \ a 2 + a X G 2 + G \ a 2 *)- Here we intro- 
duced the Pauli matrices o~ x ..z-,i referring to each qubit 
subspace, the raising (lowering) operators o~l (c^) and 
we wrote the hamiltonian in the energy basis of each 
qubit. It is more convenient to rewrite the previous 
hamiltonian in the interaction representation, resulting 
in H'j = cxp(i(H ql + H q2 )t/h)Hi cxp(-i(H ql + H q2 )t/h). 
We obtain 



Hj = hgo[exp(i(uix + U2)t)crf cr} 

+ exp(j( — U>1 — U>2)t)(Ti <T2 

+ exp(i(wi — UJ2)t)o-i &2 



(1) 



As soon as |Ai — A2I > go, the corresponding evolution 
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operator only contains rapidly rotating terms, prohibit- 
ing any transition to take place. This is a mere conse- 
quence of energy conservation : two coupled spins can 
exchange energy only if they are on resonance. 

More elaborate coupling strategies than the fixed linear 
coupling have been proposed 0, 0, . In these the- 
oretical proposals, the coupling between qubits is medi- 
ated by a circuit containing Josephson junctions, so that 
the effective coupling constant can be tuned by varying 
an external parameter (such as, for instance, the flux 
through a SQUID loop). Nevertheless, these schemes 
also require the two qubits to have the same resonant 
frequency Ai = A2 if they are to be operated at their 
optimal biasing point. If on the other hand Ai ^ A2, 
only the so-called FLICFORQ scheme proposed recently 
by Rigeti et al. [13 to our knowledge provides a work- 
able 2-qubit gate. Application of strong microwave pulses 
at each qubit frequency A.; induces Rabi oscillations on 
each qubit at a frequency vr :i . When the condition 
vr\+i>r2 = (Ai — A 2 ) is satisfied, the two qubits are put 
on resonance and they can exchange energy. It is then 
possible to realize any two-qubit gate by combining the 
entangling pulses with single-qubit rotations. Note that 
in order to satisfy the above resonance condition, the two 
qubits should still be reasonably close in energy to avoid 
prohibitively large driving of each qubit which could po- 
tentially excite higher energy states or uncontrolled en- 
vironmental degrees of freedom. Single-qubit driving fre- 
quencies of order 250MHz have been achieved for charge- 
and flux-qubits |14ll22j . In order to implement the FLIC- 
FORQ scheme, one would thus need the resonance fre- 
quencies of the two qubits to differ by at most 500MHz, 
which seems within reach for charge-qubits but not for 
flux-qubits. 

While in the scheme proposed by Rigeti et al. quan- 
tum gates are realized with a fixed coupling constant g, 
our scheme relies on the possibility to modulate g by 
varying a control parameter A. This gives us the pos- 
sibility of realizing two-qubit operations with arbitrary 
fixed qubit frequencies, which is particularly attractive 
for flux-qubits. We first assume that we dispose of a 
"black box" circuit realizing this task, as shown in figure 
\I]p, actual implementation will be discussed later. Our 
parametric coupling scheme consists in modulating A at 
a frequency ui/2ir close to Ai — A2 or Ai + A 2 . Suppos- 
ing A(t) = A + SXcosujt leads to g(t) = g + Sgcos(ujt), 
with go = ,g(Ao) and Sg = (dg/d\)S\. Then, if to is close 
to the difference in qubit frequencies ui — u>i — u>2 + S12 
while \Si2\ << \l<Ji — lj^I, & few terms in the hamiltonian 
^will rotate slowly. Keeping only these terms, we obtain 

H'j = h(Sg/2)(exp(i5i 2 t)a^ a\ 
+ exp(— ib\it)aX ovT) 

Modulating the coupling constant g allows therefore to 
compensate for the rapid rotation of the coupling terms 
which used to forbid transitions in the fixed coupling 
case, and opens the possibility to realize any two-qubit 



gate. For instance, in order to perform a SWAP gate, one 
would choose (5 = and apply a microwave pulse for a 
duration At = l/(2Sg). One could implement the "anti- 
Jaynes-Cummings" hamiltonian H'j — h(8g /2)(ct 1 " + 
a^a^) as well by applying the microwave pulse at a fre- 
quency Ll> — LOl + U>2 ■ 

We note that a recent article also proposed to apply 
microwave pulses at the difference or sum frequency of 
two inductivel y co upled flux-qubits in order to generate 
entanglement [l5|. However the proposed approach is 
ineffective if the two flux-qubits are biased at their flux- 
insensitive point. In our proposal, modulating the cou- 
pling constant between the two qubits instead of applying 
the flux pulses directly through the qubit loops overcomes 
this limitation. 

A specific attractiveness of our scheme is that the effec- 
tive coupling constant driving the quantum gates Sg/2 is 
directly proportional to the amplitude of the microwave 
driving. Therefore the coupling constant can be made 
in principle arbitrarily large by driving the modulation 
strong enough, although in practice each circuit will im- 
pose a maximum amplitude modulation and modulation 
speed which have to be respected. This is very similar 
to the situation encountered in ion-trapping experiments 
[l6|. and in strong contrast with the situation encoun- 
tered in cavity quantum electrodynamics experiments. In 
the latter, the vacuum Rabi frequency, fixed by the dipolc 
matrix element and the vacuum electric field [l2l[l3> se ^ s 
a maximum speed to any two-qubit gate mediated by the 
cavity. We also note that the coupling can be switched 
ON and OFF at nanosecond speed, as fast as the switch- 
ing of Rabi pulses for single-qubit operations. 

The hamiltonian J5| is only approximate because it 
simply omits the fixed coupling term go<J x i(J X 2- In or- 
der to go beyond this approximation, we separate the 
time-independent and the time-dependent parts of the 
coupling hamiltonian by writing Hi = Hiq + Hj(t), 
where H I0 = goa x icr X 2, Hi(t) = Sg cosLO M wta x ia X 2, and 
lomw is the frequency of the modulation. We diagonalize 
H q i + H q 2 + Hiq and rewrite Hi(t) in the energy basis 
of the coupled system (dressed states basis). We go to 
second order of the perturbation theory and use the ro- 
tating wave approximation Hiq ~ go(&i &2 + cr i~ cr 2^)- A 
complete treatment is also possible but would only make 
the equations more complex without modifying our con- 
clusions. In this approximation, denoting the coupled 
eigenstates by \i, j') = 0, 1) and their energy by E^, 
we obtain that 
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The new energy states are slightly energy-shifted com- 
pared to the uncoupled ones. However, it is remarkable 
that this energy shift does not depend on the state of the 
other qubit, since for instance E' 10 — E' 0Q = E' lt — E' 01 = 
h(A 1 +g%/{A 1 -A 2 )) = h(Ai + Su). This implies in par- 
ticular that no conditional phase shift occurs that would 
lead to the creation of entanglement. We now write 



Hi{t) = Sgcos(u} M wt)o- x ia x2 



Sg cos(wmw*){[1 



t go 

V Ai-A 2 



) 2 ]|01')(10'| +h.C 



+ 2 S ^_(|10')(10'|-I01')(01'|) 
+ |00')(11'| + |11')(00'|} 

(4) 

Writing Hi(t) in the interaction representation with 
respect to the dressed basis as we did earlier in the un- 
coupled basis shows that the presence of the coupling 
go modifies our previous analysis as follows : 1) If one 
wants to drive the 1 01') — > 1 10') transition, one needs to 
modulate g at the frequency (E[ — E Q1 )/h = A 2 — Ax + 
2po/(Aa — Ax). 2) The effective coupling constant is then 
reduced by a factor 1 — (g /(Ai — A 2 )) 2 . 3) Besides the 
off-diagonal coupling term, the time-dependent hamilto- 
nian contains a longitudinal component modulated at the 
frequency ujmw- Similar terms appear in the hamilto- 
nian of single charge- or flux-qubits driven away from 
their symmetry point and have little effect on the system 
dynamics. Driving of the 1 00') — > |11') would be done in 
the same way as discussed earlier. We conclude that our 
scheme provides a workable two-qubit gate in the dressed 
state basis for any value of the fixed coupling go. How- 
ever the detection process is simpler to interpret if the 
two-qubit energy states of Hio are little entangled, that 
is if ,g << |Ai - A 2 |. 

One might be worried that the circuit used to modulate 
the coupling constant opens additional decoherence chan- 
nels. We therefore need to estimate the dephasing and 
relaxation rates. Dephasing by 1/f noise seems the most 
important issue. In particular the need to use Joseph- 
son junction circuits to make the coupling tunable might 
be a drawback since it is well-known that they suffer 
from 1/f noise. We suppose that A = A + n(t), where 
n(t) is a fluctuating variable with a 1// power spectrum. 
From equation J3J we see that the coupling hamiltonian 
go~x\o~ X 2 gives rise to a frequency shift bv of qubit 1 res- 
onance frequency, and — 8v of qubit 2. Noise in the cou- 
pling constant thus translates into noise in the qubit en- 
ergy splittings. We now compute the sensitivity coef- 
ficients D\ z = |(00'|aff/<9A|00') - (10'|d#/3A|10')| = 
2nd(8v) / dX of each qubit to noise in A, using the frame- 
work and the notations established in [H| . We obtain 



, go d 9 
J Ai - A 2 dX 



(Ac) 



(5) 



Therefore, if go « |A 2 — Ax|, it is possible to have 
a large value of dg/dX allowing rapid operation of the 
two-qubit gate, while keeping D\. z small. In particular, 



if go = 0, the qubit is only quadratically sensitive to 
noise in A since D\ z = 0. This situation is a transpo- 
sition of the optimal point concept 0] to the two-qubit 
case. Therefore our scheme provides protection against 
1/f noise arising from the junctions in the coupling cir- 
cuit, whereas if the qubits were tuned into resonance with 
DC pulses as proposed in [Tol[llllT2] | 1/f noise would be 
more harmful. 

Given the form of the interaction hamiltonian, it is 
clear that quantum noise in the variable A can only induce 
transitions in which both qubit states are flipped at the 
same time, i. e. |0x,0 2 ) — > |lx, I2) or |lx,0 2 ) — > |0 1; 1 2 ). 
The damping rates for each transition can be evaluated 
with the Fermi golden rule similar to the single qubit 
case, and will depend on the nature of the impedance 
implementing the coupling circuit. We discuss two dif- 
ferent cases, one where A shows a flat power spectrum 
and one where it is peaked. If the coupling circuit acts 
as a resistor R thermalized at a temperature T, the re- 
laxation rate is 



r x =4ir 3 (dg/dX) 2 \h(Lu res )[ 



, YlLOr 



2tt 



[coth{- 



2kT 



1]R 



(6) 

where h(u>) = (dX/dV)(u>) is a transfer function relat- 
ing A to the voltage across the coupling circuit V. The 
frequency ui res refers to 27r(Ai + A 2 ) or 27r(Ai — A 2 ), 
depending on the transition considered. This rate can 
always be made small enough by designing the circuit in 
order to reduce the transfer function |/i(o;)|, in a similar 
way as the excitation circuits for single-qubit operations. 
In the second case we may use a harmonic oscillator with 
an eigenfrequency uj c , weakly damped at a rate k by cou- 
pling to a bath at temperature T . Now the variable A is 
an operator representing the degree of freedom of the ID 
oscillator. Therefore we can write that A = Ao(a+a'). In 
the laboratory frame, the total hamiltonian now writes 
H = H q i + H q2 + H c + Hj, where H c = hu c (a'a) and 
Hi = [go + Sg (a + a))\a xl a x2 with 5g = {dg/dX)X . 
Going in the interaction representation with respect to 
H = H q i + H q2 + H c , it can be seen that the coupling- 
contains terms rotating at lu± ± uj 2 ± ui c . Thus as soon 
as the eigenfrequency of the coupling circuit is close to 
Ai± A 2 , the qubit eigenstates will be mixed with the har- 
monic oscillator states. This is certainly not a desirable 
situation if one wishes to "simply" entangle two qubits. 
Even if uj c ^ Ax ± A 2 , there will be a remaining damping 
of the qubits via the coupling circuit yielding a relaxation 
time of the order [(to c — uj res )/5go] 2 k -1 , where again u> res 
refers to 27r(Ai + A 2 ) or 27r(Ai — A 2 ) depending on the 
transition considered. In addition, fluctuations of the 
photon number induced for instance by thermal fluctu- 
ations may cause dephasing |20J if huj c is comparable to 
kT. Given all these considerations, it seems desirable 
that the frequency uj c be as high as possible, and far 
away from the qubit frequencies. We note that this sim- 
ple analysis would actually be valid for any control or 
measurement channel to which the qubit is connected, 
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FIG. 2: a) Circuit proposed in [l2T | to implement a tunable 
coupling between two flux-qubits. The two qubits are di- 
rectly coupled by a mutual inductance M qq , and also via the 
dynamical inductance of a DC-SQUID which depends on the 
bias current at fixed flux bias. The total coupling constant 
g is shown in b) for the same parameters as were considered 
in [T^ as a function of s = Ii/2Iq. The dashed line indicates 
g — 0. Inset : (dg/ds) as a function of s. 



and therefore does not constitute a specific drawback of 
our scheme. 

We will now discuss the physical implementation of the 
above ideas. Simple circuits based on Josephson junc- 
tions, and thus on the same technology as the qubits 
themselves, allow to modulate the coupling constant at 
GHz frequency [Hill. 

To be more specific in our dis- 
cussion, we will focus in particular on the scheme dis- 
cussed in [l2T ]. and show that the very circuit analyzed 
by the authors (shown in figure^) can be used to imple- 
ment our parametric coupling scheme. Two flux-qubits of 
persistent currents I q .i and energy gaps A, (i = 1, 2) are 
inductively coupled by a mutual inductance M qq . They 
are also inductively coupled to a DC-SQUID with a mu- 
tual inductance M qs . The SQUID loop (of inductance 
L) is threaded by a flux $5 , and bears a circulating cur- 
rent J. The critical current of its junctions is denoted 
Iq. Writing the hamiltonian in the qubit energy eigen- 
states at the flux-insensitive point, equation (2) in [l2j 
now writes H = -(h/2)(A 1 a z i + A 2 a z2 ) + hgo xl a x2 , 
where g = (M qq \I ql I q2 \ + M 2 qs \I ql I q2 \Re(d J/d$ s ) h )/h. 
In figure we plot the coupling constant g as a func- 
tion of the dimensionless parameter s = 1^/21$ for the 
same parameters as in 12] : Iq = 0.48/iA, L = 200pH, 
I q i = I q 2 = 0.46uA, M qq = 0.25 P H, M qs = 33pH, 
$ s = 0.45$o- We see that g strongly depends on s. 
In particular g(so) = for a specific value sq. On the 
other hand the derivative dg/ds is finite (for instance, 
dg/ds(so) = 1GHz) as can be shown in the inset of figure 
I2t>- Biasing the system at so protects it against 1/ / flux- 
noise in the SQUID loop and noise in the bias current. 
At GHz frequencies, the noise power spectrum of s is 
ohmic due to the bias current line dissipative impedance, 
and has a resonance due to the plasma frequency of the 
SQUID junctions. This resonance is in the 40GHz range 
for typical parameters and should not affect the coupled 
system dynamics. 

As an example, we now describe how we would gen- 
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FIG. 3: Calculated evolution of the density matrix under the 
application of an entangling microwave pulse at the frequency 
Ai — A2I in the SQUID bias current, with Ai = 5GHz, 
A 2 = 7GHz, g(t) = go + Sg cos(2tt(A 2 - Ai)t) and Sg = 
100MHz. For the black curve, go — ; for the grey curve, 
go = 200MHz. 



erate a maximally entangled state with two flux qubits 
biased at their flux-noise insensitive points, assuming 
Ai = 5GHz and A 2 = 1GHz. We fix the bias current 
in the SQUID to I\, = 2sqIo and start with the ground 
state |0i , O2} ■ We first apply a ir pulse to qubit 1 thus 
preparing state |li, O2). Then we apply a pulse at a fre- 
quency A2 — Ai = 2GHz in the SQUID bias current of 
amplitude 5s = 0.015. This results in an effective cou- 
pling of strength 8g/2 = (dg/ds) (5s/ 2) = 50MHz. A 
pulse of duration St = 5ns suffices then to generate the 
state (|0i, 2 ) + |li, l 2 ))/V2. We stress that thanks to the 
large value of the derivative dg/ds, even a small modula- 
tion of the bias current of Sib = 2IqSs = l5nA is enough 
to ensure such rapid gate operation. We performed a 
calculation of the evolution of the whole density matrix 
under the complete interaction hamiltonian g(t)a x \a x2 
with the parameters just mentioned. We initialized the 
two qubits in the 1 1 1 , 02) state at t — ; at t — 10ns 
an entangling pulse g(t) — Sg cos 27r(Ai — A 2 )t and last- 
ing 20ns was simulated. The result is shown as a black 
curve in figure We plot the diagonal elements of the 
total density matrix. As expected, poo,oo(i) = Pii,n = 0> 
and 010,10 = 1 - An, 01 = (cos(27r((5g/2)i)) 2 . We did 
another calculation for the same qubit parameters but 
assuming a fixed coupling go — 200AJ Hz. Following 
the analysis presented above, we initialized the system 
in the dressed state 1 10') and simulated the application 
of a microwave pulse g(t) = 5 g cosu Mwt at a frequency 
Wjfw = 2.04GHz taking into account the energy shift 
of the dressed states. The evolution of the density ma- 
trix elements (grey curve in figure shows that despite 
the finite value of go, the two qubits become maximally 
entangled as previously. The evolution is not simply si- 
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FIG. 4: (Color online) Flux qubit parametrically coupled to 
an LC oscillator via a DC SQUID, a) Electrical scheme : 
the qubit (blue loop) is inductively coupled to a DC SQUID 
shunted by a capacitor and thus forming a LC oscillator, b) 
Dependence of the coupling constant g\ as a function of the 
bias current [2p|l . At the current II the coupling constant 
vanishes, c) Derivative dgi/dlt as a function of It stays 
nearly at a constant value on the current range considered. 



nusoidal because we plot the density matrix coefficients 
in the uncoupled state basis. Note also the slightly slower 
evolution compared to the go = case, consistent with 
our analysis. This shows that the scheme should actually 
work for a wide range of experimental parameters. 

It is straightforward to extend the scheme discussed 
above to the case of a qubit coupled to a harmonic os- 
cillator of widely different frequency. As an example 
we consider the circuit studied in [20L l2ll |22| which 
is shown in figure A flux qubit is coupled to the 
plasma mode of its DC SQUID shunted by an on-chip 
capacitor C s h (resonance frequency v p ) via the SQUID 
circulating current J. As discussed in |2j}, the cou- 
pling between the two systems can be written Hi = 
[gi(I b )(a + a" 1 ) + g 2 (Ib)(a + a^) 2 ]a x . We evaluated gi(h) 
for the following parameters : $5 = 0.45<I>o, h = lfJ-A, 



qubit-SQUID mutual inductance M = lOpH , qubit per- 
sistent current I p = 240nA, A = 5.5GHz, v p = 9GHz 
as shown in figure At J& = 1% = 0, the coupling con- 
stant gi vanishes. It has been shown in [2(| that when 
biased at lb = 1^ and at its flux-insensitive point, the 
flux-qubit could reach remarkably long spin-echo times 
(up to 4/is). On the other hand, the derivative of g\ is 
shown in figure to be nearly constant with a value 
dgi/dlb — ~4:GHz/fiA. Therefore, inducing a modula- 
tion of the SQUID bias current Si cos(2n(vp — A)t) with 
amplitude 5i = 50nA would be enough to reach an effec- 
tive coupling constant of 100MHz. The state of the qubit 
and of the oscillator are thus swapped in 5ns for reason- 
able circuit parameters. This process is very similar to 
the sideband resonances which have been predicted [24[ 
and observed j2^. However, in order to use these side- 
band resonances for quantum information processing, the 
quality factor of the harmonic oscillator must be as large 
as possible, contrary to the experiments described in |23 | 
where Q ~ 100. This can be achieved by superconducting 
distributed resonators for which quality factors in the 10 6 
range have been observed |23| . Employing this harmonic 
oscillator as a bus allows the extension of the scheme to 
an arbitrary number of qubits, each of them coupled to 
the bus via a SQUID-based parametric coupling scheme. 

In conclusion, we have presented a scheme to entangle 
two quantum systems of different fixed frequencies cou- 
pled by a a x a x interaction. By modulating the coupling 
constant at the sum (difference) of their resonance fre- 
quencies, we recover a Jaynes (anti-Jaynes) -Cummings 
interaction hamiltonian. It also yields an intrinsic pro- 
tection against 1/f noise in the coupling circuit. Our pro- 
posal is well suited for qubits based on Josephson junc- 
tions, since they readily allow tunable coupling constants 
to be implemented. The idea can be extended to the in- 
teraction between a qubit and a harmonic oscillator and 
could provide the basis for a scalable architecture for a 
quantum computer based on qubits, all biased at their 
optimal points. 
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